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Synopsis 

Discrete particle simulations by Accelerated Stokesian Dynamics (ASD) and a microstruc¬ 
tural theory are applied to study structure and the viscosity of hard-sphere Brownian sus¬ 
pensions in active microrheology (MR). The work considers moderate to dense suspensions, 
from near to far from equilibrium conditions. The microscopic theory explicitly considers 
many-body hydrodynamic interactions in active MR, and is compared with the results of 
the ASD simulations, which include detailed near and far held hydrodynamic interactions. 
We consider probe and bath particles which are spherical and of the same radius a. Two 
conditions of moving the probe sphere are considered: these apply constant force (CF) and 
constant velocity (CV), which approximately model magnetic bead and optical tweezer ex¬ 
periments, respectively. The structure is quantihed using the probability distribution of 
colloidal particles around the probe, Tb|p(r) = ng{r), giving the probability of hnding a bath 
particle centered at a vector position r relative to a moving probe particle instantaneously 
centered at the origin; n is the bath particle number density, and is related to the suspension 
solid volume fraction, (j), hy n = Sb/dvra^. The pair distribution function for the bath parti¬ 
cles relative to the probe, g{r), is computed as a solution to the pair Smoluchowski equation 
(SE) for 0.2 < (j) < 0.50, and a range of Peclet numbers, describing the ratio of external force 
on the probe to thermal forces and dehned Re/ = (fc^T) and Peu = (kbT) 

for CF and CV conditions, respectively. Results of simulation and theory demonstrate that 
a wake zone depleted of bath particles behind the moving probe forms at large Peclet num¬ 
bers, while a boundary-layer accumulation develops upstream. The wake length saturates 
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at Pe/ ^ 1 for CF while it continuously grows in CV. This contrast in behavior is related 
to the dispersion in the motion of the probe under CF conditions, while CV motion has no 
dispersion; the dispersion is a direct result of many-body hydrodynamic interactions. This 
effect is incorporated in the theory as a force-induced hydrodynamic diffusion flux in the 
pair SE. We also demonstrate that, despite this difference of structure in the two methods 
of moving the probe, the probability distribution of particles near the probe is primarily set 
by the Peclet number, for both CF and CV conditions, in agreement with dilute theories; 
as a consequence, similar values for apparent viscosity are found for the CF and CV condi¬ 
tions. Using the microscopic theory, the structural anisotropy and Brownian viscosity near 
equilibrium are shown to be quantitatively similar in both CF and CV motions, which is in 
contrast with the dilute theory which predict larger distortions and Brownian viscosities in 
CV, by a factor of two relative to CF microrheology. This difference relative to dilute theory 
arises due to the determining role of many-body interactions associated with the underlying 
equilibrium structure in the semi-dilute to concentrated regime. 


1 Introduction 

The field of microrheology (MR), as an alternative approach for measuring rheology of com¬ 
plex fluids at microscopic length scales, has undergone tremendous progress over the past two 
decades [MacKintosh &: Schmidt, 1999; Waigh, 2005; Squires & Mason, 2010; Wirtz, 2009]. 
The motion of a probe particle (or particles), subjected to internal stimuli in passive MR or 
external stimuli in active MR, is measured and used to characterize the material properties in 
the vicinity of the probe(s). One appealing aspect of MR is that it allows analysis of the lo¬ 
cal microscopic structure of the material studied, facilitating investigation of microstructural 
heterogeneity [Valentine et al, 2001; Dasgupta & Weitz, 2005]. Another advantage of MR 
over bulk rheology is the significantly smaller volumes required, and its non-intrusive nature 
makes it ideal for rheological study of living systems [MacKintosh &: Schmidt, 1999; Wirtz, 
2009]. Other advantages include the large accessible frequency range in viscoelastic response 
[Willenbacher & Oelschlaeger, 2007], and direct measurement of interactions [Verma et ai, 
1998]. 

In a pioneering work. Mason &: Weitz [Mason &: Weitz, 1995] demonstrated that the linear 
viscoelastic response of the material can be extracted by monitoring the thermal fluctuations 
of a micron-sized probe particle embedded in the complex fluid using the generalized Stokes- 
Einstein relation (GSER). The theory developed by Mason &: Weitz relates viscoelastic 
moduli to statistical measures of the probe motions such as mean square displacement or 
velocity auto-correlation [Squires & Mason, 2010]. Two major assumptions of the theory are 
i) the fluid must be near equilibrium and ii) the motion of the fluid about the probe should 
be inertialess, i.e. described by Stokes equations. A number of conditions can violate these 
two assumptions [Squires & Mason, 2010; Squires, 2008; Depuit & Squires, 2012a,6]. We 
focus on strongly nonlinear conditions, as these play a central role in the rheology of complex 
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fluids. In many applications, complex fluids are driven far from equilibrium by externally 
applied stresses, resulting in nonlinear responses including shear thinning and thickening, as 
well as normal stress differences. Far from equilibrium, the GSER cannot be applied. Active 
microrheology is an alternative method to investigate the nonlinear response by externally 
moving the probe through the fluid using, for example, magnetic forces [Habdaas et al, 2004] 
or optical tweezers [Sriram et ai, 2010]. Important questions addressed here are: To what 
degree can the results of active MR be compared against bulk shear rheology? What are 
the main differences between the two methods, and do those differences offer a new insight 
about the material? How is the behavior measured dependent upon the specific method 
of generating the probe particle motion? Specifically how do the measured structure and 
viscosity differ if the probe’s velocity is controlled externally (CV) relative to the case where 
the probe is moved with a controlled external force (CF)? 

The hard, or near-hard, sphere colloidal (Brownian particle) suspension is ideal for in¬ 
vestigating these fundamental aspects of active MR, as a theoretical foundation of some 
completeness is already available for these systems, while the simulational tool of Stokesian 
Dynamics provides access to the structure and stresses to guide theoretical development for 
MR. Previous theoretical studies of MR for colloidal suspensions can be divided into two 
groups. The first group of works, by Brady and coworkers, are microscopic theories devel¬ 
oped for dilute systems (cj) ^ 1) with interactions limited to pair level. For example. Squires 
& Brady [2005] studied dilute colloidal dispersions with only hard-sphere interactions (ne¬ 
glecting hydrodynamic interaction), and computed the distribution of bath particles around 
the probe as a function of the ratio of the external force to thermal forces (Peclet number, 
Pe). They observed that as the external force increases the distribution of colloidal particles 
around the probe becomes more anisotropic, with large accumulation of particles in front 
and a depleted zone or wake behind the probe. This work was extended to consider the 
effect of pair hydrodynamic interactions in dilute suspension [Khair &: Brady, 2006], the role 
of probe particle’s shape [Khair &: Brady, 2008], relation between stress and force-induced 
diffusivity of the probe [Zia &: Brady, 2012], and the case of having an oscillating external 
force/velocity [Swan et al, 2014]. The predictions of the dilute theories are in qualitative 
agreement with available experimental work [Sriram et al, 2010] and Brownian Dynamics 
simulations [Carpen & Brady, 2005]. The viscosity computed from the microstructure also 
shows qualitatively similar behavior as a function of Peclet number to shear viscosity al¬ 
though there are clear quantitative differences. A second group of theoretical studies has 
focused on dense suspensions near the jamming transition. These theories ignore hydrody¬ 
namic interactions between the probe and colloidal particles and use mode-coupling theory 
[Gazuz et al, 2009; Voigtmann & Fuchs, 2013] to predict the microstructure and from this 
extract the MR of colloidal suspensions. While these methods are powerful for studying very 
dense systems, due to the complete absence of hydrodynamic interactions, the hydrodynamic 
stress is not represented in these models. This limits the ability to accurately predict non¬ 
linear behaviors, such as shear thickening in shear flow and the analogous force-thickening 
in active MR. There has also been a recent simulational study of active micrheology of soft- 
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particle systems near jamming transition [Mohan et al, 2014], The structure and rheology in 
this system are determined by elastohydrodynamic contact forces, which are not considered 
in our study. 

In this study, we focus on the range of parameter space that has not been previously 
studied for MR, namely concentrated (up to ^ = 0.50) near-hard sphere colloidal suspen¬ 
sions and Peclet numbers ranging from near- (Pe <C 1) to far-from-equilibrium [Pe 3> 1) 
conditions. We use Accelerated Stokesian Dynamics (ASD) simulations [Banchio & Brady, 
2003] to study structure (pair distribution function, g{r)) and rheology (apparent viscosity) 
by modeling the active microrheology experiments. This is, to our knowledge, the hrst re¬ 
port on the discrete particle simulation of colloidal suspensions in active MR that includes 
detailed near- and far-held hydrodynamic interactions. Two conditions of regulating the mo¬ 
tion of the probe are considered: the constant force (CF) approach, controlling the external 
force, on the probe which mimics the experimental approach of pulling a magnetic bead 
probe [Habdaas et al, 2004; Choi et al, 2011]; and the constant velocity (CV) approach of 
setting the velocity, of the probe, as is to a close approximation done by optical tweez¬ 
ers [Wilson et al, 2009; Sriram et al, 2010]. The Peclet number for CF and CV conditions 
are respectively dehned as Pej = {F''^^a)/khT and Peu = {67rr]LP'^^a'^)/kbT, where a is the 
radius of the probe, ry is the viscosity of the huid medium, and is the thermal energy. 
The cases of CF and CV microrheology may be regarded, respectively, as the analogs of 
bulk rheometry in controlled stress and controlled rate of strain modes, which typically yield 
similar results. In active MR it is not certain that the two methods are equivalent. 

We use sampling from ASD simulations to show that the structural anisotropies induced 
by the motion of the probe in CF and CV are distinctly different at Pcf^u 1: In CV the 
depleted zone of bath particles behind the probe monotonically extends with increasing the 
Peclet number (or external velocity). By contrast when moving the probe with constant force, 
the length of depleted zone (and other structural features away from the probe) saturate with 
Pef. We demonstrate that this saturation is a result of the dispersive motion of the probe in 
CF conditions, while in CV the probe’s trajectory is a straight line. We show that the source 
of the dispersion in CF is the hydrodynamic interactions of the probe with its neighboring 
bath particles. This is the hrst study that takes into account the effect of hydrodynamic 
dispersion on the structure and viscosity of colloidal suspensions in active microrheology; 
the previous theories ignore this effect by either assuming the suspension is dilute and the 
interactions are limited to pair particles [Squires & Brady, 2005; Khair & Brady, 2006; Swan 
&: Zia, 2013; Swan et al, 2014], or the HI is entirely ignored [Gazuz et al, 2009; Voigtmann 
& Fuchs, 2013]. 

We complement the simulation results with a microscopic theory developed by modifying 
our theoretical framework for shear hows [Nazockdast &: Morris, 2012&,a, 2013] for active 
microrheology. We consider both conditions of CF and CV. The theory is based on com¬ 
putation of the probe-bath particle pair correlation, described through a pair distribution 
function, g{v), determined as a solution to the appropriate form of the pair Smoluchowski 
equation (PSE). This theory captures through probabilistic means of incorporating many- 
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body effects the interactions that lead to dense-packed structural features, e.g. next-nearest 
neighbor rings. The dispersion of the probe’s trajectory in CF conditions is incorporated 
in our theory by adding a diffusive flux term in the PSE, which points to the other unique 
aspect of this theory. The force-induced diffusion coefficient is modeled by extending the 
equilibrium relationship between diffusion and particle pressure to non-equilibrium condi¬ 
tions. The predictions of the theory are in good agreement with the ASD results over the 
range of Pe and (f) studied, although discrepancies will be noted. 

Our simulation results and predictions of theory show that the boundary layer structure 
of the density of the particles formed very near the probe at Pe 1 scales similarly with Pe 
for CF and CV (unlike the structure away from the probe). We specifically show that the 
pair distribution function is determined by the net external force on the probe, irrespective 
of the volume fraction and method of moving the probe. This finding is in agreement with 
the predictions of the dilute theories. Due to this similarity of near-contact structure for CF 
and CV conditions, similar values of g were computed for apparent viscosity at Pe 2> 1 from 
theory and simulation for the two forcing conditions. 

The near-equilibrium structure and rheology is only studied theoretically, due to the poor 
sampling statistics of ASD simulations at Pe < 1. The predicted structural distortion at 
Pe <C 1 for CF and CV are similar in magnitude and form, resulting in similar values for 
their respective viscosities. This is not in agreement with predictions of dilute theories which 
predict structural distortions and Brownian viscosities twice as large for CV conditions. We 
show that this difference arises because of the near-contact excluded volume interactions 
which integrate to an effective mean force from the bath particles in PSE; this effect is 
completely absent in the dilute theories where the pair distribution function is near unity 
everywhere. These many-body interactions result in a significantly faster decay of distorted 
structure away from the probe in the present theory as compared to the predictions of the 
dilute theories. 


2 Simulation of active microrheology 

The Accelerated Stokesian Dynamic (ASD) simulation method is applied to compute the 
dynamics of N Brownian hard-sphere particles interacting hydrodynamically. The sampled 
positions of the spheres are then used to compute the pair distribution function, g{v), and 
make direct comparison with predictions of the theory. The details of the conventional 
Stokesian Dynamics (SD) and the more computationally efficient implementation ASD are 
given elsewhere [Sierou & Brady, 2001; Banchio & Brady, 2003]. A brief description of the 
method and its implementation for active microrheology are discussed here. 

Stokesian Dynamics computes the dynamics of particles immersed in an incompressible 
Newtonian fluid with negligible inertial effects, i.e. at vanishing Reynolds number. The 
hydrodynamic interaction between the particles is governed by Stokes equations and the 
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motion of the particles is determined through Langevin equation 

m . + F^^^ (1) 

at 

where in a system composed of N particles, m is the generalized mass/moment-of inertia 
6N X 6N matrix, U is the generalized translational/rotational velocity vector of particles of 
size 6N, and F^, F^, F'^ are generalized hydrodynamic. Brownian, interparticle force/torque 
vectors of size 6N, respectively. The Langevin equation is then integrated in time steps large 
enough that the momentum of the particles is relaxed (m • dU /dt 0) At ^ = m/Q'Krja, 

but small enough that the conhguration of the particles is not substantially changed due to 
Brownian diffusion or external flows; the result is the over-damped Langevin equation for 
Brownian particles. This equation is again integrated in time with At 3> r™, to obtain the 
evolution equation for the positions and alignments of the particles [Banchio & Brady, 2003]. 

In active microrheology, a force is also applied externally to drag the probe through 
the bath particles which is represented with F®^h Note that is zero for all particles 
except the probe. In the absence of externally imposed flows, we have F^ = • U 

for hydrodynamic and interparticle forces where is the grand resistance tensor of size 
QN X QN which relates force to velocity of the particles and only depends on the positions of 
the particles. For particles near contact, the resistance matrix is computed using a pairwise 
additive lubrication resistance formulation and when particles are well-separated, a pairwise 
additive mobility formulation is adopted, with the mobility matrix inverted to obtain the 
resistance matrix. The Brownian force is random with zero time average, and strength given 
according to the fluctuation-dissipation theorem by (F^(0)F'®(t)) = 2kbT'R^'^6(t). Once the 
resistance matrix is known from the conhguration of particles, the velocity of the particles 
including the probe can be computed as a solution to (1) . Afterwards, (1) is integrated in 
time to update the positions of particles. Periodic boundary conditions are considered in all 
three directions. 

The implementation for constant force on the probe (CF) is straight forward; the forces, 
including applied only on the probe, are explicitly known at each time step, with F^ 
drawn from a random distribution. The velocities are calculated as a solution to the set 
of linear equations dehned by (1). Under precise CV conditions, F-^ (again, on the probe 
alone) as well as and TJ-® are unknowns, which presents difficulties that were avoided 
by numerically replicating what is done in laser tweezer experiments [Sriram et al, 2010]. 
Instead of pulling the probe with a truly constant velocity we simulated the reverse 
case of hxing the probe while the bath particles are all in the external flow of = — U®"^*. 
A spring force proportional to deviation from its initial position was applied to the probe, 
F™*(t) = —A'(xi(t) — Xi(0)) where K is the spring constant. As the simulation begins the 
probe is displaced from its initial position in the direction of the external flow until the 
spring force, in average, balances the force exerted by the bath particles. Due to fluctuations 
in microstructure, the probe position fluctuates around its average value, which is slightly 
displaced in the flow direction, |Ax|/a ~ /K. This slight fluctuation may prove 
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Figure 1: The trajectories of the probe particle in CV condition from simulations with different 
initial conditions at Peu = 10 and cj = 0.40, projected to x — z plane, where z is the direction 
of the probe motion. 


important relative to a precisely fixed velocity, but is actually a better physical mimic of 
experimental conditions. As an example of the behavior, in hgure 1 we illustrate trajectories 
of the probe under CV conditions at Peu = ° = 10 and (f = 0.40. These are 

obtained from simulations with different initial conditions. The trajectories are projected 
to the X — z plane where z is the direction of the probe’s motion; note that x and y are 
equivalent directions owing to axisymmetry about the probe’s mean motion. As pointed 
out, the probe undergoes a mean displacement in the direction of the flow of a magnitude 
(|x|/a) ~ 0.6, while the displacement in other directions is zero in average. Viscosity is 
computed as yu = — (F®’’’*)/(67rl7®’’’*a) where = — K ((xi) — Xi(0)) is the ensemble 

average of the external force, taken after the displacement of the probe reaches a statistically 
steady behavior. The spring constant was chosen proportional to Peu at Peu > 1 such that 
probe remains within a radius of 0.2a from its initial position for the large majority of the 
time. 

The number of configurations needed to achieve a statistically meaningful structure in 
active microrheology is significantly larger than in shear flow as the structure is measured 
only around the probe. Thus,we took statistics from 30 to 60 simulations over an extended 
period of motion for each Pe and (f. The structural distortions from isotropy induced by 
pulling the probe through the suspension near equilibrium {Pe 1) are proportional to 
Peclet, and this was not found to be feasible to access by simulation with our available 
resources. To illustrate the problem, our simulation results from 60 simulation runs at 
Pef = \, (j) = 0.40 yield < 0.07 where is the maximum value of g{v). 

For smaller Pe/, the magnitude of anisotropy will be smaller, and we find that 

the number of simulations performed in this study is not large enough to obtain a statistically 
reliable measure of g — gegi^) at Pe/ < 0.5. Thus we present ASD results for Pef >1 and 
compare with the theoretical results, while for study Pe/ < 1 only the theory is used (see 
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§ 6 . 1 ). 

The far-from-equilibrium condition, i.e. Pe 3> 1, is characterized by the formation of 
a wake region depleted of bath particles behind the probe. This wake can extend to many 
particle diameters, resulting in a complication for the simulation technique applied here, 
as the probe may interact with the wake formed by its periodic image. To avoid this, 
simulations with N = 2048 particles were performed for these conditions. Nevertheless, 
under CV conditions we hnd the wake to be especially pronounced, and in fact it extends 
over the entire length of the simulation box (over ten particle diameters) at Peu > 25. 

3 Simulated microstructure and trajectories 

In this section, we present a limited set of results, considering only the simulated probe 
trajectories and the microstructure at Pe > 1. These are obtained from ASD simulation 
implemented as described in the prior section. The microrheology viscosity results deter¬ 
mined from simulation (as described in §5) will be presented in §6 and 7, where they will be 
compared with the theory. The theory is formulated in the following §4. The results of this 
section serve to emphasize the key difference between CF and CV conditions. In the latter 
where its velocity is fixed, the probe travels in an essentially straight line through the sus¬ 
pension, with only the slight fluctuations allowed by the spring force. In contrast to this, for 
a fixed force, the probe changes its direction and its speed in response to the heterogeneity 
of the local structure of the suspension bath. We begin by showing that this difference in 
the trajectory of the probe results in distinct large-scale structural differences when Pe 3> 1. 
Here, we recall that Pe denotes the ratio of the external force applied on the probe to ther¬ 
mal forces and will be specified to Pe/ = F'''^^a/{kbT) and Peu = /{ki,T) for CF 

and CV conditions, respectively, where needed for clarity. Here, rj is the viscosity of the 
suspending fluid. 

The microstructure is described by the distribution of number density of bath particles 
around the probe, P(r) = ng{r) where n = N/V is the average number density and g{r)dr 
is thus the normalized likelihood of observing a bath particle between r and r -|- dr; note 
that this normalization implies g' —>■ 1 at large r = |r| where the influence of the probe on 
the bath particle distribution vanishes. 

Figures 2(a) to 2(c) illustrate ^(r) from sampling the simulated configurations of the bath 
particles around a probe (whose center serves as the origin) at ^ = 40 and Peu = 10, 25 
and 100. At Peu = 10 a boundary layer structure is formed with a zone of large ^'(r) near 
contact in the direction of motion and a zone depleted of bath particles (wake zone) behind 
the probe. The length of the wake zone increases with the flow strength, and at Peu = 25 
the wake spans the entire simulation box. Consequently, the probe interacts with its periodic 
image which results in a depleted streak in front of the probe at Peu = 25 and Peu = 100. 

The sampled structures for CF condition at (/) = 0.40 and Pe/ = 100, 400, and 1000 are 
shown in figures 2(d) to 2(f); recall Pe/ = . A boundary layer structure, similar to that 




(a) Peu = 10, (^ = 0.40, ASD 


(b) Peu = 25, (l> = 0.40, ASD 


(c) Peu = 100, 4 , = 0.40, ASD 



ASD 



(g) P^f = 100 


(h) Pef = 400 


(i) Pef = 1000 


Figure 2: (a-f) The computed pair micro structure, ^(r), 0 / bath particles around the probe 
obtained from ASD simulations; (a-c) constant velocity mode of motion at <j) = 0.40 and Peu = 
10, 25, and 100, respectively; (d-f) constant force mode at f = 0.40 and Pcf = 10, 100, and 
1000, respectively. The maximum value in each color bar indicates the scale of the near-contact 
pair probability. The lower limit of the color bars representing 0(1) features of g{r) are similar 
for all Peclet numbers to facilitate comparison. The panels 2(g) to 2(i) illustrate simulated 
trajectories of the probe, projected onto the x-z plane worth z the direction of motion, under 
CF conditions at cf = 0.40 and Pef = 100, 400 and 1000, respectively. A dispersion around the 
average trajectory of the probe is observed, with similar dispersion for Peg = 400 and 1000. 
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found in CV conditions, is observed at Pe/ = 100. The near contact values of ^'(r) in front 
of the probe (shown as the maximum value on the color bars) monotonically increase from 
Pcf = 100 to 1000. The length of the wake zone increases from Ptf = 100 to Pcf = 400. 
Aside from the increase of contact values with Peclet number, the structure (including the 
length of the wake zone) remains unchanged from Pej = 400 to 1000. This saturation of the 
structure far from contact, with the most noticeable feature being the wake length, marks the 
essential difference we have observed between CF and CV microrheology (MR) ai Pe ^ 1. 

In hgures 2(g) to 2(i) we present the trajectories of the probe particle at Pef = 100, 
400 and 1000, corresponding to the microstructures in hgures 2(d) to 2(f). The physical 
origin of the contrast between CV and CF structure at Pe 1 is seen in these trajectories, 
which clearly show a hydrodynamic dispersion around the average motion of the probe. The 
observed dispersion at Pef ^ 1 is a consequence of hydrodynamic interactions (HI), and 
does not have thermal origin. This becomes evident by noting that the variation of the 
dispersion radius with traveling distance remains almost unchanged as the conditions are 
changed from Pef = 400 to Pef = 1000, i.e. the dispersion is proportional to the distance 
traveled, regardless of the time taken to travel the distance. The wake zone is hlled 

in by particles owing to the hydrodynamic dispersion, reaching uniform number density 
(g 1) at sufficient distance from the probe: this distance saturates at large Pe because 
the dispersion itself saturates. The observations show that the length of the wake zone 
at Pe 3> 1 is roughly the axial distance traveled in order for the probe motion to show 
a signihcant probability of lateral displacement by its own radius through the dispersion 
process; the length of the wake zone eventually saturates as Pe —oo and for (j) = 0.40, the 
wake extent is approximately 10a. For CV conditions, there is only limited hydrodynamic 
dispersion of the bath particles and not of the probe, so the wake hlls in primarily by thermal 
diffusion of the bath particles, which is much slower than the translation of the probe at 
large Peu, leading to the very long wakes at hxed probe velocity. 

Because the dispersion of the probe’s trajectory is induced by many-body HI, it increases 
as (j) increases, resulting in reduction of the length of the wake. We demonstrate this by 
showing g{r) and probe trajectories at Pe/ = 100 and different volume fractions in hgure 3. 

A quantitative measure of the dispersion around the average motion can be obtained by 
computing the mean-square displacement (MSD) in the plane perpendicular to the direction 
of motion. In hgure 4(a) we present the MSD as a function of the distance travelled, denoted 
probe at (/) = 0.40 and Pe/ = 100, 400 and 1000; this is taken for convenience 
along one coordinate direction, y, but the dynamics is axisymmetric around the probe’s 
motion (along z), and the dispersion in x is equivalent. Note that for distances traveled 
smaller than the particle size, the extent of dispersion is larger at the smaller values of 
Pe/. For > 5, however, we see convergence of the displacements at all Pe/ to a 

single curve. The observed behavior can be explained by noting that the dispersion at small 
-g inhuencod by Brownian motion, which at a hxed velocity is stronger at smaller Pe/. 
The force-induced dispersion is a result of the probe interacting with a sequence of different 
conhgurations of neighboring particles as it moves through the suspension, and when the 
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(a) Pef = 100, 0 = 0.30, ASD (b) Pef = 100, </> = 0.45, ASD (c) Pef = 100, </> = 0.50, ASD 
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(d) Pef = 100, </> = 0.30, ASD (e) Pef = 100, <j) = 0.45, ASD (f) Pef = 100, ij> = 0.50, ASD 


Figure 3: (a-c) The sampled g{r) from ASD simulation under CF conditions at Pcf = 100, 
for (f) = 0.30, 0.45 and 0.50. The wake zone is reduced in size with increasing (f. (d-f) The 
trajectories of the probe for individual simulations for the same conditions as figures (a-c), 
respectively. The dispersion around the average motion increases with volume fraction resulting 
in a reduction of the wake zone. 



Figure 4: The mean squared displacement (MSD) of the probe particle motion, with distances 
scaled by the particle radius a. (a) MSD in the directions orthogonal to the direction of motion 
at (f) = 0.40 and Pcf = 100, 400, and 1000; (b) MSD at Pcf = 100 and (f = 0.30, 0.40, 0.45 
and 0.50. 
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probe has interacted with many bath particles > 10a), the force-induced dispersion 

dominates the Brownian diffusion such that the extent of dispersion shown as MSD vs 
becomes independent of Pe. We plot the MSD for different volume fractions in figure 4(b), 
showing the increase of dispersion at larger (j). 

We note that the force-induced dispersion has not been included in previous theories for 
active microrheology of colloidal dispersions. The previous theories have either neglected 
many-body interactions by assuming the suspension is dilute [Squires & Brady, 2005; Khair 
&: Brady, 2006; Swan & Zia, 2013], or only consider steric interactions and simply discard 
hydrodynamic interactions [Gazuz et ai, 2009; Gnann et ai, 2011; Voigtmann &: Fuchs, 2013]. 
In the next section we present our microscopic theory which incorporates the force-induced 
dispersion as well as the average velocity induced by many-body particle and hydrodynamic 
interactions. These features are found to be necessary to allow the theory developed here to 
reasonably predict the microstructure and viscosity up to cj) = 0.40 and Pe 2> 1. 


4 Microstructural theory formulation 

Here we develop the theoretical description, which follows the approach used in our theory 
developed for the microstructure of sheared suspensions [Nazockdast &: Morris, 20126, 2013]. 
We describe first the case of a constant velocity probe, and then the constant force case. 
Results of the theory will be presented in §6 and §7, with comparison to ASD results at 
Pe > 1. 

4.1 Constant velocity probe 

We consider N identical spherical particles dispersed in an incompressible Newtonian fluid. 
One of these particles, labeled 1, is defined to be the probe particle that is pulled through the 
bath of particles with a constant velocity (GV), Ui. The ultimate goal is to compute the force 
exerted on the probe by the bath particles and fluid. The motion of the probe particle through 
the suspension changes the spatial distribution of force- and torque-free bath particles around 
it. The resistance to motion of the probe depends on this distribution, and consequently the 
force and suspension microstructure are coupled. The Ai-body configurational probability 
distribution function, P^(xi, X 2 ,..., Xjv), satisfies the conservation equation 

OPn . 

^ Vq, • Jq, — 0, 

0.-1 

where is the probability flux associated with particle a. When pulling the probe with 
a fixed velocity, the probability flux associated with the probe is purely convective and is 
expressed as ji = UiP/v- The remaining particles move freely in response to Brownian, 
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interparticle and hydrodynamic forces, with the flux for a ^ 1 expressed as 

//3=N \ /3=N 

j. = ^ . Fn Fa, - ^ . V^Fa. (3) 

\/3=l ) 13=1 

Here is the hydrodynamic velocity of particle a induced by motion of the probe particle 
through its hydrodynamic interaction with bath particles, is the mobility of particle 

a due to a force on particle and is the Brownian diffusion tensor. The 

velocities of particle a induced by interparticle and Brownian forces applied on the rest of 
the particles are given, respectively, by 

j3=N (3=N 

Uj = ■ Fj and Uj = - ^ D.a ■ V^Fa/Fa. 

13=1 13=1 

To allow progress, we reduce Pn to the distribution of number density of bath particles 
around the probe, P{v) = ng{r) where n = N/V is the average number density and ^'(r) is 
the pair distribution function for the bath-probe pair. To reduce to an equation for g{v), (2) 
is integrated over all possible conhgurations of — 2 particles keeping particle 1 and one 
bath particle (which we refer to as particle 2) hxed (meaning of known position, but mobile). 
This yields 

|+V(j2-j.>,=0. (4) 

The average of a function a4(r) at a hxed separation is denoted (H) 2 (r), with the averaging 
performed over all possible conhgurations of the other N — 2 bath particles. We make a 
mean-held approximation, assuming that the particle pair of interest is immersed in a held 
that contains the average effect of the other particles. As a consequence, we neglect the 
convective hux associated with huctuations of force and mobility: 


(M.F),= (M),.{F),. 


(5) 


Applying the mean-held approximation to (4), the average velocities of particles 1 and 2 are 
related to the average hydrodynamic forces on each particle as follows: 


(Ff) 


H\ 1 
2 


where 


.(Ff) 

2- 

_R21 

Rii 

Ri2 

r 

_R21 

R 22 . 



R-ii R-121 r ~Ui 

-(U2)i 


-1 


( 6 ) 


( 7 ) 
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Note that the formulation is written for hydrodynamic forces applied from the fluid to the 
particles, resulting in a negative sign on velocities in (6). Since the probe velocity is assumed 
hxed, (Uf )2 = 0 and the Brownian force on particle 2 from the bath particles is 

(Ff )2 = -hTVg{v)/g{v). (8) 


Recall that the bath particles are force- and torque-free on timescales longer than the mo¬ 
mentum relaxation time, (F 2)2 = (^ 2)2 + ^ 2)2 + ^ 2)2 ~ incorporating (6) and (8) 

to this condition for the bath particles, (U 2)2 is written as 


(U2)2 = (R22)-' • 


• Ui + {F^), - hTVg{v)/g{r) 


(9) 


Taking ( j2 — ji )2 = ((U 2)2 — Ui) (?(r), (4) in steady state is modihed to the pair SE for CV 
conditions. 


Vr 


C ■ Ui^(r) + M • {F^),g{r) - hTM ■ Vrg{r) 


= 0 , 


( 10 ) 


where C = — (R22) ^ • R21 — I and M = (R22) Quantities are made dimensionless as 
follows: 


r-a, (M„^)2 -Mo, Ui - (FO2 - 

where Mq = -is the single particle mobility. The pair SE in dimensionless form is 


V- 


C • Ui^(r) + • {F^),giv) - V^(r 


= 0 , 


( 11 ) 


and we recall that 

Uexta ^ Q^^uext^2 

~ hTMo ~ ■ 

The boundary conditions for (11) are 

jr = 0 at r = 2, 

g ^ 1 as r —)■ 00 . 


( 12 ) 


(13a) 

(13b) 


The hrst boundary condition imposes zero radial flux at contact due to excluded volume. In 
order to solve (11), the conditional averages in the equation including average pair mobility 
and interparticle force should be formulated in terms of r and ^(r). The formulations are 
given in §4.3. 
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4.2 Constant force probe 


We now consider the case where the probe particle is pulled through the bath with a constant 
force (CF), with the force taken along the z axis. The dispersive flux induced by 

velocity fluctuations of the probe and bath particles is the fundamental difference of CF 
and CV active microrheology. In a CF experiment, unlike when the probe is dragged with a 
constant velocity, the probe can change direction and speed. As a result, the probe trajectory 
is not always in the direction of the applied force, resulting in the noted dispersion. If the 

d 

probe is pulled at Pe/ = ^ 1, thermal fluctuations become negligible with respect 

kbT 

to the force-induced fluctuations in setting the long-time dispersion; thus it is essential to 
include dispersion in the probe motion when formulating the structure problem for active 
mirorheology with CF at Pef ^ 1. 

The CF formulation follows the CV formulation until the relation (4). The probability 
fluxes are calculated based on a mobility formulation for CF conditions. The hydrodynamic 
velocities of the probe and particle 2 are expressed as 


■(Uf); 




"jpext" 

im'if), (Mg,f)2. 


. 0, . 


( 14 ) 


which gives 


u" = {V^), - (Uf)^ = - ({M"'-), - (M"/),) ■ F“‘. (15) 

Note that the bath particles are force-free, i.e. = 0 and F®^* in (14) is applied externally; 

as a result, the hydrodynamic force applied on particles from the fluid to balance the external 
force is —F®^*, and in (14), unlike (6), there is no negative sign applied to the velocity 
components. The relative velocity associated with interparticle forces can be computed 
using a relation similar to (14), given by 

u'’ = {U,f)j - (Uf), = ({M22 >j - {Mi2>j) • (FP2 + ((M2,>2 - {M„)2) ■ {Ff)2, (16) 


where the (7 F superscript of the mobility tensor M is dropped for convenience. 

We model the flux associated with the flow-induced dispersion as a diffusion. The relative 
flux due to velocity fluctuations is therefore written as 

j? - jf = - (D22 - Di 2 + Dn - D21) ■ V^(r), ( 17 ) 

where 


D 22 — ^ 6 ^( 1 ^ 122 ) 2 ) (18a) 

Di 2 = ^6^(1x412)2, (18b) 

Dll = ^6^(1x411)2 -|-D{i, ( 18 c) 

D 21 = fc6T(M2i)2 + D 21 , (18d) 
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where the D superscript on j in (17) is indicative of the diffusion ffux. The force-induced 
dispersion only appears on D 21 and Du, since only interparticle and thermal forces (and 
no external force) act on particle 2 (D^ and 0(2 are zero); recall that Dq ,/3 is the diffusion 
of particle a induced by the force on particle jS. The diffusion induced by Brownian forces 
scales as k^TM^ where kbT/a is the thermal force, and Mq is the single particle mobility. 


When Pcf = 


kbT 


^ 1, the diffusion of the probe is expected to scale with instead of 


kbT/a, i.e. D-^ MoF^^^a. Based on the given relations, the CF pair SE reduces to 


Vr 


(U^(r) + - (D^(r) + D-^(r)) • Vr^(r) 


0 , 


(19) 


with 

D"" = fc,T((Mn )2 - (M 2 i )2 + (M 22)2 - (Mi 2 ) 2 ), (20a) 

D'/'(r) = D{i(r)-D4(r), (20b) 

where D'®(r) and D'^(r) are relative diffusivities induced by Brownian and external forces. 
Quantities are made dimensionless using the scalings 

r-a, - fcbTMo/a, - MoF®^*a. 


The dimensionless form of the pair SE is 


V- 


(^U^(r)-FPe^^U^(r))^(r) - ^D^(r) + D-^(r)) ■ V^(r) 


0 , 


( 21 ) 


where Pcf = 


1 ^' 


The boundary conditions are the same as the CV case, given by (13). 


The average mobility and interparticle force for the probe and particle 2 should be formulated 
in terms of r and g{v). In addition D'^ also needs to be defined. These issues are addressed 
immediately below in §4.3 and §4.4. 


4.3 Average mobility and interparticle forces 

The formulation of average mobility and interparticles forces is as given in our formulation for 
simple-shear flow which is explained in detail in earlier works [Nazockdast h Morris, 20126, 
2013]. Here we only give the final expressions. There is, however, one key difference between 
conditional averages in active microrheology and shear flow. In the shear flow, the motion is 
induced by an external straining flow, a uniform field acting on each individual particle; by 
contrast, in microrheology only the probe particle is subjected to an external force/velocity. 
As a result, when the Smoluchowski equation is integrated to pair particle level in shear 
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flow, particles 1 and 2 are statistically identical, resulting in certain important symmetries: 
in shear flow (Mii )2 = (M 22)25 (^ 12)2 = (]V[ 2 i) 2 , and (F 2)2 = —(Fi) 2 - On the other hand, 
in the pair SE for active microrheology (whether for CF or CV), the distribution of bath 
particles around the probe is intrinsically different from that around particle 2 since the force 
or velocity is only imposed on the probe. A clear example would be the case where particle 
2 and the probe are well separated. In this configuration, particle 2 is decorrelated from the 
probe; hence particle 2 and the bath particles around it are in the equilibrium distribution, 
fl'(r 32 ) = 5 '®®(?" 32 ) where r 32 is the separation vector of particle 2 and a third bath particle. 
On the other hand, as it is shown in previous works and will be shown here, the distribution 
of bath particles around the probe is anisotropic with more particles accumulating in the 
direction of probe motion and a depleted zone of particles behind the probe. As a result the 
symmetry properties in shear flow do not hold here, and in general 

(Mn)2 / (M22)2, ^ -(Ff)2. 

This effect of the lack of symmetry between the probe and bath particle of the pair becomes 
negligible for dilute suspensions, since the interactions are limited to the pair level (particle 
1 and 2) while for concentrated suspensions it should be included in the formulation. This 
is taken into account through modification of the triplet distribution function, g' 3 (r, ri 3 ), as 
described below. 

Interparticle forces are assumed to be pair-wise additive and the conditional average 
forces are formulated as 


(Fi)2^(r) = F2i(r)^(r) +n 


F3i(r3i)^3(r,r3i)*3i, 


( 22 ) 


where the forces F 21 and F 31 are exerted on particle 1 by particles 2 and 3, respectively. The 
hard-sphere force is modeled here through a steep repulsion written as 


Ff 


-ihT/a) 


rexp(—re) 

----—r 

1 — exp(—re) ’ 


(23) 


where e = r/a — 2 is the dimensionless gap between the particle surfaces. The parameter 
r determines the steepness of the repulsive force; larger r results in a steeper decay of the 
repulsion, and we have used r = 400 for all the calculations here. It is necessary to express 
the triplet distribution, ^ 3 , in terms of the pair distribution function to obtain a relationship 
which is only a function of ^'(r) and r. This is done here through a modified Kirkwood 
superposition closure approximation Rice Sz Lekner [1965] given by 


93(r,r3i) = 9(r)9(r3i)9(r32)exp(r*(r,r3i,</i)). 


(24) 


This modifies the Kirkwood superposition through the exponential term which is tabulated 
based on separation distance (r 3 i here) and various values of (j) [Rice & Lekner, 1965]. To 
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take into account the fact that conhgurations of bath particles around particle 2 and the 
probe are generally different from one another this closure is modihed to 


^3(r,r3i,r32) = ^(r)^(r 3 i) (A/(r 32 ) + (1 - A)/«(r 32 )) exp(-r*), (25) 


where A is a measure of connectivity of particles 2 and 3 with the probe (particle 1). This 
is modeled as 


-'(r, r3i,r32) 


max ( U2 


IU3I) 


|Ui 



( 26 ) 


where Ui, U 2 , and U 3 are velocities of particles 1, 2, and 3 at each triplet conhguration. 
This expression gives a correct description of the limits where particles are touching or well- 
separated. For example, if both 2 and 3 are well separated from the probe particle, A —0 and 
the distribution of bath particles around particle 2 is the equilibrium distribution, g^^{r^ 2 )- 
Also when all three particles are connected in the direction of the applied force A —1 
since three particles move as a single body along the force direction. Unlike the shear flow, 
fl'(r 23 ) 7 ^ <?(r 32 )- To correctly account for this, in (26) the separation vector of particles 2 
and 3 is measured from the particle closer to 1 to the more distant particle. For example, 
in the case of having particle 3 between particles 1 and 2, ^*(r 32 ) = 5 '(r 32 ); when particle 
2 is between 1 and 3, g*{Y^ 2 ) = fl'(—r 32 )- In other words, between particles 2 and 3, the 
one closer to the probe carries the information on the motion of the probe to the particle 
at a greater distance. We note that this is a simple phenomenological model. Formulations 
of the average resistance and mobility are identical to those given in our previous paper 
[Nazockdast & Morris, 2013] and are not reproduced. 


4.4 Force-induced diffusion 

f f 

Variations of D — D 23 with separation are modeled in a similar fashion to that used in our 
theory for structure and rheology in shear flow [Nazockdast &: Morris, 20126, 2013]: 

= Df{v) = D{ ( 6 ;(r)ff + n{r) (I - ff) ), (27) 

where is the isotropic portion of long-time self-diffusion of the probe, = tr(D{)/3 
which scales as at Pe/ S> 1. The functions Q{r) and 'H(r) determine the variations 

of force-induced relative diffusivity, D'f, along the line of centers and the plane perpendicular 
to the pair radial vector, i.e. the local [6 — (p) plane, respectively. As was done for sheared 
suspensions, we model the relative dispersion based on the behavior of in the limits of 
a pair that is well-separated, r/2a 2 > 1, and very near contact, r —)■ 2a. When the pair are 
well-separated, the particles become decorrelated and the induced dispersion on particle 2 
due to motion of the probe vanishes, D 23 (r) = 0, and D{j^(r) —)► D{ where D{ is the long 
time average self-diffusivity of the probe. The relative dispersion, here a relative diffusivity, 
is proportional to relative velocity fluctuations, i.e. D-^(r) cx r(U'(r)U'(r)) 2 , where U' is 
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the deviation of the local relative velocity from the average value and r is the characteristic 
time scale over which relative motion becomes decorrelated. When the probe and particle 2 
are in near-contact conhgurations, the relative radial velocity and its fluctuations approach 
zero due to the excluded volume effect, U oc r — 2a and U' oc r — 2a. Thus the relative 
diffusivity in the radial direction is expected to approach zero as (r — 2a)The excluded 
volume restriction is imposed along the line of centers and the fluctuations of relative relative 
velocity in the angular direction remain on the order of the probe velocity. As a result, at 
all separations we take D-^ oc MoF®^*a for the 6 and (p components. 

Under equilibrium conditions, the collective diffusivity is related to osmotic pressure by 

D®® = — [Dhont, 1996] where subscript c refers to collective diffusion and mobility. 

Since mobility and pressure are still well-dehned outside equilibrium, one may extend the 

as 

equilibrium relationship to non-equilibrium systems by analogy and write D = M • where 

on 

S is the stress induced by the particle phase. Microstructural theories for predicting single 
particle and collective diffusion of sheared suspension also hnd a general correspondence 
between the shear induced stress and diffusion [Brady & Morris, 1997; Leshansky & Brady, 
2005; Leshansky et ai, 2008]. [Zia & Brady, 2012] extended this idea to active microrheology 
to relate probe diffusivity to the force-induced stress on the probe. They considered a phase 
composed of a suspension of probe particles where the suspension is so dilute that the probes 
have no interaction among themselves. The momentum balance is written for this phase and 
an exact relationship between the force-induced stress and diffusion for a dilute suspension 
is derived. Their analysis, however, neglects hydrodynamic interactions between the probe 
and the bath particles which makes a direct extension of their theory to the problem under 
study difficult. Instead of their detailed formulation, we use the analogy in the simplest form 
to compute the limiting value of D7(r) when the probe and particle 2 are well-separated, 
Df and write 


= Df = M 


dE 

dnj 


(28) 


where Ua is the number density of the probe phase ( <C 1) and S is the average stress on 
the probe phase per unit volume. Note that since the probe phase is considered very dilute 
(no interactions between the probes), the collective diffusivity and mobility become identical 
to their single particle values in (28). 

In the next step, we take the enhanced diffusion to be dominated by hydrodynamic 
velocity fluctuations (at Pe/ 3> 1 where the effect is pronounced) and hence only consider 
the hydrodynamic contribution to the stress in (28). The reasonable assumption that the 
hydrodynamic stress is dominated by near-contact lubrication stress allows S to be written 
in terms of the moment of external force in the radial direction: 

l>r* 

= UaU / . r) ffg'(r)(ir. (29) 

J2a 
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The integration of stress is carried out inside the boundary layer of ^'(r) from contact to 
r = r*, where r* is taken to be the radius at which the increase in angularly averaged pair 
correlation function, g{r), from the equilibrium distribution drops to 10% of the contact 
values i.e. g{r*) = 0.1 [g{2a) — g^'^{2a)] + g^‘^[r*). This is similar to the dehnition used by 
Sierou & Brady [2002] as well as our earlier work on shear flow [Nazockdast & Morris, 2013]. 
Note that D/ has no dependence on Ua after taking the derivative with respect to Ua in (28). 
From (28) and (29) the force-induced diffusivity is found to scale roughly as D{ ~ U'^^^(j)ag{2) 
where g{2) is the contact value of pair distribution function in the boundary layer portion 
of the contact surface. 

For simplicity we assume the variations of relative diffusion with (j) take the same form as 
the far-held value, D(, for all pair separations. As a result Q and T-L are only functions of r. 
The near-contact relative radial diffusion caused by the force on the probe varies as (r — 2)^ 
(with lengths scaled by a) near contact and approaches the force-induced self diffusion when 
the probe and particle 2 are well separated. Based on these limits Q is modeled as 

Q{r) = 1 — exp [—m(r/a — 2)^] , (30) 

where m is a constant for all conditions of (f) and Pef\ this is the only htting parameter in 
our theory. Setting rh = 0.01 resulted in a good agreement between theory and simulation. 
We did not attempt to systematically vary to hnd the optimal value of fh. 

Since there is little experimental or simulation data on active microrheology to provide 
meaningful averages of the probe and bath dynamics, we use the insight from shear how 
Accelerated Stokesian Dynamics simulations for modeling the variations of angular diffusiv¬ 
ity. The results of simulation sampling at Pe 3> 1 show the magnitude of shear-induced 
angular relative velocity huctuations remains hnite and is reduced to approximately 0.25 of 
the maximum far-held value [Nazockdast & Morris, 2013]. We assume these limits hold and 
model 71 as 

'^(^) = i 

The conditional averages appearing in (11) and (19) are expressed in terms of integrals 
of g{r). The resulting pair Smoluchowski equations for both CF and CV conditions are 
nonlinear integro-differential equations which we solve using a hnite element iterative scheme 
with the boundary conditions given in (13). The details of the solution technique are given 
in Nazockdast & Morris [20126]. There is, however, one difference in computation domain 
of the active microrheology and shear how. In shear how, the structure changes in all three 
directions in space while active microrheology is an axisymmetric problem with symmetry 
about the direction of external force on, or velocity of, the probe. 
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5 Viscosity 

In this section, we develop the relationship between ^'(r) 
calculated using the Stokes-Einstein relationship, rjmicro 
F — pext -g goal compute as a function of (?(r), implying the need 

for the average probe mobility, (Mii) 2 . For CV, the problem is inverted: the velocity is 
hxed and the external force should be calculated which requires computation of the average 
resistance to motion of the probe. 


and apparent viscosity which is 
F 

= —. In a CF experiment, 
Gvrt/a 


5.1 Constant velocity viscosity, rju 

The motion of the probe at any separation from particle 2 due to Brownian motion is the 
sum of the forces directly applied from the fluid, (Ff )2 = kbTVg/g and the hydrodynamic 
interaction owing to the Brownian motion of particle 2 given by — R 12 • where Uf = 
—kbT(R. 22 )~^ ■ '^g/g is the Brownian velocity of particle 2. Combining these relationships, 
the total Brownian force on the probe simplihes to 


= n kbT 




Ri2 • (R: 


■22 J 


\-il 


[Vg/g) ■ zg{r)dr = n kbT / (V • C) • zg{r)dr, 


(32a) 


where subscript z implies the component of the force along the direction of probe motion 
and C = —R 22 ^ • R 21 — I. Integration by parts and the divergence theorem were used to recast 
the second integral form from the hrst expression on the right hand side for in (32a). 
similar formulation is used to obtain the net interparticle force on the probe resulting in 


A 


F[ = -nkbT / [C ■ Ffi(r)] • z^(r)dr 


(32b) 


Finally the hydrodynamic force applied from the fluid on the probe is simply the product of 
the average resistance of the probe in the direction of motion and its velocity: 




(32c) 


The summation of the forces in (32) gives the force applied on the probe from the fluid and 
bath particles. Thus Fy* = — [F^ + F^ + F^^ needs to be externally applied to the probe 
to maintain its constant velocity of The integrals in (32) are easily computed once g{v) 
is known from the theory. The micro-viscosity is calculated as 


riu_ 

rj 


= Vu = 


- (rf + Ff + F») 

67rry[/®^*a 


(33) 


We only consider hard-sphere interparticle interactions. Thus, F^ is only non-zero at con¬ 
tact. The hard-sphere force is modeled using a steep repulsive radial force, (23). Using 
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simple scaling arguments, it is straightforward to show that Ff oc Hence the inter¬ 

particle force is neglected here in comparison with Brownian forces. This assumption does 
not generally hold for particles interacting through long-ranged soft potentials and the con¬ 
tribution from interparticle interaction needs to be computed Nazockdast & Morris [2012a]. 
The hydrodynamic dimensionless viscosity for hard-sphere interaction is thus 

fi§ = {Rll)/Ro, (34) 

where Rq = Qn'qa is the single particle resistance. The Brownian micro-viscosity is fj^ = 
flu - Vu- 

5.2 Constant force viscosity, r]f 

The CF micro-viscosity is dehned as 

P'ext 

where is the average velocity of the probe particle in the direction of the applied 

force. The average probe velocity is stated in terms of contributions from hydrodynamic, 
interparticle and Brownian interactions, i.e. 


(35) 


+ u^. 


(36) 


with the contributions related to g{v) through the following relations [Khair h Brady, 2006]: 


= (Mil) • Ff*, 

V^ = nj (Mil - Mi 2)2 • Fg{r)g{r)dr, 
= -n / V • (Mil - Mi 2 ) 2 ( 7 (r) dr. 


(37a) 

(37b) 

(37c) 


The coefficients appearing in (37a - 37c) are known through solving for ^(r). For near 
hard-sphere dispersions, the contribution of bath particles to viscosity is taken to be a 
summation of contributions from pure fluid. Brownian, and hydrodynamic interactions. The 
hydrodynamic contribution to viscosity is dehned as 


rj {^n )' 


(38) 


where (Mff) is the average probe mobility in the force direction. The Brownian viscosity is 
obtained by subtracting hydrodynamic viscosity from the total viscosity: f}^ = f}f — fjf- 
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6 Predictions of the theory 

In presenting the structure at weak Pe, or near-equilibrium conditions, the theoretical results 
presented in §6.1 stand alone because statistically meaningful results were not accessible in a 
reasonable time by simulation. The results of theory and simulation at Pe > 1 are compared 
in §6.2. 

6.1 Near-equilibrium, Pe <C 1 

Near equilibrium, the distortion of microstructure from its equilibrium form can be formally 
expressed as a regular perturbation expansion in terms of Pe as 


g^{^) = geq{r) + Peugi (r) -F 0{Pel), (39a) 

/(r) = geq{r) + Pe/^((r) + 0{Pe}), (39b) 

with 

gi (r) = fi {r)U^^* 2 . ■ f (39c) 

gi (r) = f{ {r)F^^^z ■ f, (39d) 


where subscripts U and / present the constant velocity and force conditions. When Pe —)► 0, 
0{Pe^) contributions to the structure become negligible, and the zero-force/velocity viscosity 
is determined by and fj. To calculate f( and we compute g{r) for (f) < 0.40 at Pe/ and 
Peu = 0.01, for which contributions of 0{Pe^) and higher order terms can be neglected. Note 
that (j) = 0.40 is the largest volume fraction for which the modihed Kirkwood superposition 
gives a convergent solution for ^'(r) [Rice & Lekner, 1965]. The expressions (39c) and (39d) 
are then htted with good agreement to the predicted microstructures at Pe/,;/ = 0.01 to 
give fi{r) and f( (r). To study the effect of CF and CV conditions, the structures should 
only be compared against each other when they both correspond to the same external force 
or velocity. The choice of = (Rn) • results in This constraint is 

satished by rescaling the Peclet number for CF as Pe(;7) = Pe/(Mff)“^ which is analogous 
to comparing fi{r) against f (where Mff = M^^/Mq. 

Figure 5(a) shows the variations of /^(2) and /j(2)(Mff)“^ with (j), showing that the 
values increase with volume fraction and the distortions are slightly stronger in CV. The 
difference is stronger at higher (j). Having roughly similar structural distortion in the CF 
and CV scenarios is not in agreement with the predictions of Squires & Brady [2005] for 
dilute Brownian suspensions with no hydrodynamic interactions (HI), in which a distortion 
twice as large for constant velocity compared to constant force was predicted. Recall that 
both probe and the bath particles diffuse due to Brownian fluctuations in CF while in 
CV only the bath particle does so. This gives rise to a factor of two enhancement to the 
pair diffusion in CF compared with CV when HI are ignored causing the distortion of pair 
microstructure in CV to be twice as large as CF conditions. At finite volume fractions, the 
presence of many-body interactions change this simple analysis. In the equilibrium limit. 
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Figure 5: (a) The contact values ofO{Pe) structural deformation functions for CF and CV, 
{Mii)~^fj{2) and /^(2) in (39a, b), vs 4>; (b) the variations of fj with radial separation, r, at 
different volume fractions compared against the predictions of the dilute theory (equation (4-9) 
in Khair & Brady [2006]). 


these interactions introduce a potential of mean force, whose effects include the secondary 
peaks in pair correlation function, and nonlinear variations of material properties such as 
the viscosity, the osmotic pressure, and the diffusion coefficient with volume fraction. Below 
we give a brief explanation of these effects and demonstrate how they alter the structure at 
Pe <C 1 with respect to the predictions of dilute theories. 

In figure 5(b) we compare the radial variations of f( at different volume fractions against 
the predictions of the dilute theory (equation (4.9) in Khair & Brady [2006]). The main 
observation is that the decay of structural distortion away from contact is significantly faster 
in the predictions of the full theory than the dilute theory, even at the lowest volume fraction 
studied, (j) = 0.20. The structural buildup at contact from the present concentrated theory 
drops rapidly and oscillates about zero for r > 2.5, whereas in the dilute theory the distortion 
is monotonically decreasing and is still 1/4 of the value at contact when r = 4. The secondary 
peaks in f( are generated by propagation of the anisotropic pair correlation buildup near 
contact through interactions with the bath particles [Nazockdast & Morris, 20126]; a similar 
propagation of correlation in equilibrium results in the secondary peaks in geq{r). The sharp 
decay of structural anisotropy compared to the dilute limit, and the presence of secondary 
peaks in fj both show that the near contact interactions determine the structure at volume 
fractions as low as (f = 0.20. 

The extent of distortion of equilibrium structure at small Pe is proportional to the ratio of 
convective to diffusive fluxes. In figure 6(a), we compare the dimensionless relative diffusion 
coefficients from CF, 2 (Mii — M21), and CV, M = R 22 ^, at ^ = 0.33. At large separations, as 
expected, the relative diffusion of CF conditions is twice as large as the CV case while in both 
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Figure 6 : (a) The variations of the dimensionless relative diffusivities with pair separation 
distance, r, under CF and CV conditions and (f = 0.33. At r ^ 1, the CF diffusion is twice 
the value of the CV diffusion while both diffusivities approach zero at contact; (b) the results 
of figure 6(a) shown for 2 < r < 2.5. The isolated pair relative diffusion coefficients from the 
dilute theory are also shown for comparison. 


cases, the values asymptote to zero at contact. In figure 6 (b) the variations near contact, 
in 2 < r < 2.5, are considered. We have also included the variations of relative diffusion 
coefficients of an isolated pair in the dilute theories for comparison. For both the dilute and 
concentrated suspensions in this range of pair separations, the ratio of CF to CV relative 
diffusivities is at most 1.2 compared to 2 in Brownian suspensions with no HI. This trend was 
observed for other volume fractions as well. Since the structural distortion predominantly 
occurs in r < 2.5 as shown by figure 5(b), where the relative diffusion coefficients for CF and 
CV are quantitatively close, our theory predicts similar distortions for the two conditions at 
a given average velocity; see figure 5(a). 

For further insight, we analyze the governing equations for g\j and The equations, 
obtained by substituting (39a) and (39b) in (11) and (21), are 

=) aL + ufmf + \j{g{ - 2M^ • = 0, 


V 

V 


[(C ■ z 
-M. 


9^q + 


(40a) 

(40b) 


where Mj, = Mu — M 12 . Similar to perturbation expansion of (?(r), we have expanded the 
velocity induced by interparticle interactions, U, as 

U/T = ^TU ^ peu{’^, (40c) 

with 


K., = M ■ 

V{„ = M, ■ (Fr - Ff)i,„ 
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(40e) 











Figure 7: The ratio of 0{Pe) convective fluxes from interparticle interactions, {\5^)igeq, 
to the flux induced by the average motion of the probe, {C ■ i)geq, under CV conditions as a 
function of radial separation; see equation 4-1. 


where C = —(R 22 )~^ • R 21 ~ I- For readability we have dropped (.)2 around conditionally 
averaged quantities. The solution to the pair Smoluchowski equation in equilibrium can also 
be used to write (Ff)e 9 = —(Ff)eg = —V logg^qipr) which implies this is the force arising 
from the well-known potential of mean force in liquid-state theory; the many-body flux terms 
in (40a) and (40b) are not present in the dilute theories. One measure of the importance 
of these particle interactions can be obtained by taking the ratio of the 0{Pe) anisotropic 
convection flux induced by interparticle interactions, (U^)^'^ -|- (Uf)( 7 ^, to the convection 
flux induced by average motion of the probe, (C ■ z)geg. The predictions of the theory show 
that generally <C \\J^g^g\. In this case the ratio of two convection fluxes in radial 

direction simplihes to 

r = ■ f) (f ■ C • z)“^. (41a) 

Due to the linearity of the structural distortion with 17®^^ the radial velocity is decomposable 
to the general form of 


r = z-fJ^(r). (41b) 

We plot ^{r) for various volume fractions in hgure 7. At all the ratio is signihcant and 
for (f> > 0.30 it becomes larger than 1 for regions near the boundary; this is indicative of the 
dominating role of the underlying structure and its effect on particle interactions in setting 
the microstructure in the near-equilibrium dense suspension. The results for CF conditions 
are quantitatively close to those of hgure 7. 
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6.2 Far from equilibrium, Pe ^ 1 

For Pe > 1, simulation results are available for comparison against the theoretical predictions 
of microstructure. Figure 8 compares the predictions of ^'(r) in the plane of symmetry with 
the results obtained from sampling conhgurations from ASD simulations at Pej = 10,100, 
and 1000. At Pe/ = 10 a clear distortion of equilibrium structure is observed. The distortions 
become stronger with Pe/, leading to a boundary layer structure on the leading half of the 
probe surface. In both theory and simulations the structure away from contact saturates at 
Pe/ = 1000 as a result of force-induced dispersion being the predominate far-held relative 
diffusivity. The predictions for Pe/ = 400 and Pe/ = 1000 are quite close except for the 
maximum values of ^'(r) at contact. We do not present the pe/ = 400 condition in hgure 
8. At all Peclet numbers, the second nearest neighbor peaks are distorted similarly to g at 
contact, i.e. the values are largest in the how direction and smallest behind the probe. There 
is a general agreement between the predictions and simulation results. However, differences 
can be noted, including the fact that the form of the theory and simulational distortion 
at Pe/ = 10 behind the probe differ. The next nearest neighbor peaks are stronger in the 
predicted microstructure, indicating that the integral formulation for many-body interactions 
overpredicts the correlation between the particles. This was also observed in our previous 
study of pair dynamics in sheared suspensions [Nazockdast & Morris, 2013]. Also, in the 
ASD results the width of the farheld wake zone narrows with increase of pair separation 
distance more than is seen in the theoretical predictions. 



(a) Pef = 10 


(b) Pef = 100 


(c) Pef = 1000 


Figure 8: The variations of g{r) in the (r — z) plane passing through the probe particle 
diameter, at fixed force from simulation (ASD) and theory, for ft = 0.40. (a) Pe/ = 10, (b) 
Pef = 100, and (c) Pef = 1000.Identical color bars are chosen for theory and simulation. 


To investigate the performance of the theory at different volume fractions, we compare 
the predictions with the sampled structures from ASD for Pe/ = 100 at (f = 0.30 and 
0.50 in hgure 9. The theory correctly predicts the reduction of the length of the wake zone 
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(a) <f) = 0.30, ASD (b) <f) = 0.30, ASD 



(c) (j) — 0.30, Theory 


(d) (j) — 0.50, Theory 


Figure 9: Theory predictions and results of simulation (ASD) for g{v) under CF conditions 
at Pcf = 100; (a) ASD and c) = 0.30, (b) theory and 4> = 0.30, (c) ASD and 4> = 0.50, and 
(d) theory and cj) = 0.50. 


with increasing (f due to the enhancement of force-induced dispersion at higher particle 
loading. In the theory, the length of the wake zone is determined by the ratio of convective 
to diffusive fluxes outside of the boundary layer. When force-induced diffusion is accounted 
for, the diffusion flux at Pej 1 is determined by D'f oc MF®^*5'(2)(/) (see §4.4). In this 
limit Jd/Jc oc (j) leading to a smaller wake at higher cf, and saturation of the wake length for 
Pef > 1 . 

Figure 10 shows ^'(r) for (j) = 0.40 and Peu = 10, 25 and 100, from both theory and 
simulation. The most prominent feature seen in CV forcing, unlike CF, is that the length 
of the wake continuously increases with Peu, and eventually the wake spans the entire 
simulation box. The length of the wake at Peu = 25 is larger than Pef = 1000 which 
corresponds to approximately Pe(^u) = = 285. Recall that this rescaling of 

Peclet number corresponds to equal velocities in both CF and CV conditions. The theory 
developed here gives reasonable predictions of structure in both CF and CV scenarios and 
shows that the observed differences can be attributed to the dispersive motion of the probe for 
CF conditions, a flux mechanism absent under CV conditions. A difference between theory 
and ASD results can also be noted by comparing figure 10(a) and figure 10(d). Similar to 
CF conditions, the theory overestimates the magnitude of the next nearest neighbor peaks. 
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(a) Peu = 10, ASD 


(b) Peu = 25, ASD 


(c) Peu = 100, ASD 
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(d) Peu — 10, Theory (e) Peu — 25, Theory (f) Peu — 100, Theory 


Figure 10: The prediction and simulation results of g{r) under CV conditions at (j) = 0.40 
and varying Peu- (a) ASD and Peu = 10, (b) ASD and Peu = 25, (c) ASD and Peu = 100, 
(d), theory and Peu = 10, (e) theory and Peu = 25, and (f) theory and Peu = 100. 


6.2.1 Microstructure near the probe 

The influence of probe motion on average resistance (or viscosity) is controlled by near¬ 
contact interactions of the probe with bath particles, which is quantified by the near-contact 
form of gij). The values of pair microstructure at contact, 3'(2;(/?), from theory and simu¬ 
lation are compared in figure 11 at Pej = 1,10 and 100 and (j) = 0.40. The theory gives a 
good quantitative prediction of simulation results in the range of Pe studied. 



Figure 11: The predictions and simulation results of pair distribution function at contact as 
a function of azimuthal angel, g{2, ip), (a) Pep = 1, (b) Pep = 10, and (c) Pep = 100. 
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Next we examine the effect of volume fraction. Figure 12(a) shows g{2-, ip), for Ptf = 100 
and (j) = 0.30, 0.40, 0.45 and 0.50 from simulational sampling and predictions of the theory. 
The variations are almost independent of volume fraction, except for a slight deviation for 
(j) = 0.30. This can be understood by analyzing the ratio of convection to diffusion terms 
in pair SE. The diffusion at contact is controlled by Brownian diffusion = —fc{,(M) 2 V( 7 . 
Since mobility appears as a multiplying factor in both the convective and diffusive flux 
terms, changing the volume fraction does not change the ratio of the two fluxes, jo/jc = 
PeJ^'Vlng{r), and near contact structure remains almost unchanged while varying ^ at a 
given 




V 


(c) CV-Normalized 


Figure 12: Simulation results and theory prediction for g{2, p) at different volume fractions 
for CF conditions at Peg = 100 and CV conditions at Peu = 10. (a) CF, (b) CV, (c) the 
results of figure 12(b), multiplied by the dimensionless equilibrium mobility of the probe at the 
given volume fraction, Meqif) = Meg{4>)/Mo. The solid lines are the predictions of the theory. 
In figures 12(a) and 12(c), only the predictions at cj) = 0.40 are illustrated. Predictions in other 
volume fractions are quantitatively very similar and are not presented to aid visualization. 


In figure 12(b) we illustrate g{2,<p) ai (f = 0.20, 0.30, 0.40 and 0.45 for CV conditions. 
To avoid interaction with the periodic image in ASD simulation, Pejj = 10 was chosen for all 
volume fractions. The predictions are presented with solid lines. There is good agreement 
between theory and simulation results. Figure 12(c) presents the same results with the pair 
distribution function multiplied by equilibrium dimensionless mobility. The curves collapse 
well on a single curve, with some deviation at (j) = 0.20. Similar to CF conditions, the near 
contact pair microstructure is controlled by the force required to pull the probe. Note that 
the force-induced diffusion becomes smaller than Brownian diffusion inside the boundary 
layer; hence we expect similar scaling of near contact microstructure with the applied force 
(or average velocity) and with cf for both CF and CV conditions. 

Figure 13 shows the maximum contact value of the pair distribution function as a function 
of Peu for CV and the equivalent Peclet number for CF based on the average probe motion, 
Ijjext)^ The equilibrium contact value of g was subtracted to consider only the structural 
distortion due to external forcing. When CF values are plotted against Pegj)/2, they lie on 
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Figure 13: Simulation results for increase in maximum contact value of pair distribution 
function from equilibrium, — 5 ™“^, as a function of Peu for CV and Pe(ij^/2 for CF. The 
solid line is the best power-law fit to the data with exponent of n = 0.95. 


the same line as CV values, i.e. the near contact microstructures of CV and CF conditions are 
matched when they are compared at equal Peu and Pe(u) /2. The factor of 1/2 was multiplied 
to Pe(^u) based on the arguments given by the previous dilute theories stating that in CV 
only the bath particles can move by Brownian motion, making the relative diffusivity half 
that of the CF conditions. Although the details of hydrodynamic interactions may change 
the shifting factor by some amount, the analysis of dilute suspensions seems to give a good 
representation of microstructure at contact. Thus we hnd that the main shortcoming of the 
dilute theories dX Pei^ \ is the absence of force-induced diffusion which results in incorrect 
predictions of pair microstructure away from contact under CF conditions. 


7 Microrheology 

Figure 14 shows the predicted dimensionless hydrodynamic. Brownian, and total micro¬ 
viscosity as a function of (f for CF and CV. The Accelerated Stokesian Dynamics simulation 
values (obtained using a Green-Kubo relation) are also presented for comparison. The pre¬ 
dicted values for zero-force/velocity hydrodynamic micro-viscosity are identical for both CF 
and CV since the average mobility near equilibrium is independent of Pe and is dehned by 
(recall that the hydrodynamic viscosity is the inverse of the average probe mobility). 
Hence only CV predictions are presented. 

It is evident that the values of Brownian viscosities are signihcantly smaller than shear 
flow in all volume fractions. We note that the Brownian viscosities for CF and CV are 
quantitatively close, with CV values becoming slightly larger with increasing (f. This is 
similar to the trend observed for 0{Pe) structural deformations in hgure 5(a). Since under 
both conditions is determined as an integral of the product of structural anisotropy 
and divergence of the pair velocity in pair space (see (32a) and (37c)), similar structural 
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distortions result in similar Brownian viscosities. We again observe a contrast between our 
predictions and the results of dilute theory theories that predict r]§ « 2r]f [Squires & Brady, 
2005; Swan & Zia, 2013]. 



(a) Brownian Viscosity 



(b) Hydrodynamic Viscosity 



(c) Total Viscosity 


Figure 14: The predicted near-equilibrium CF and CV viscosities as a function of (f com¬ 
pared against Accelerated Stokesian Dynamics (ASD) simulation results for shear flow, (a) the 
Brownian viscosity, (b) the hydrodynamic viscosity, and (c) the total viscosity given by the 
summation of the Brownian and hydrodynamic viscosities. 


Figure (14(b)) shows that the predicted near-equilibrium hydrodynamic viscosity is smaller 
than the values obtained in shear flow for all volume fractions. There is a very good agree¬ 
ment between fj^ and the solid line in Figure 14(b), which presents the empirical relation¬ 
ship for the inverse of the dimensionless short-time self-diffusion coefficient at equilibrium, 
Do/Dl{(j)) = [(1 — 1.56^)(1 — 0.27^)]~\ This result was given by Lionberger & Russel [1997] 
based on fitting to experimental measurements of Pearson & Shikata [1994] on hard-sphere 
suspensions near equilibrium, with Dq = kfTMQ the isolated particle short-time self-diffusion; 
recall that fj^ = Mq/{M) = Dq/D l{(j),Pe). This indicates that the zero force/velocity mo¬ 
bility of the probe is predicted accurately for (f < 0.40. 

Figure 15 shows variations of theoretical predictions of CF and CV viscosities with the 
relevant Peclet number for a ^ = 0.40 hard-sphere suspension along with our active mi¬ 
crorheology ASD simulation results for total viscosity. The results are compared against the 
simulation results of viscosity in shear flow at different Pca^ = ■ The ASD results for 

microrheology are only presented for a limited range of Pe due to the difficulty of getting 
statistically reliable data at low Pe and the large simulation box size needed for Pe 
particularly for CV conditions. For these limited data points, ASD results are in good agree¬ 
ment with predictions for viscosity. The hydrodynamic contributions to viscosity under both 
CF and CV predictions follow the trend observed in the shear flow results and show very 
little change with Peclet number until Pe = 1. At larger Peclet numbers, the values mildly 
increase with Pe in the range studied here. The predicted values for fj^ are very close for CF 
and CV, reflecting the similar near-contact microstructure. Brownian contributions in all 
cases are reduced with increase of Pe and become negligible at Pe > 1. The total viscosity 
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Figure 15: Viscosity as a function of Pe for cf = 0.40 at (a) constant force and (b) constant 
velocity. The open symbols are the predicted values for microrheology and filled symbols are 
the results of ASD simulations in shear flow at the equivalent Peg. The stars present the total 
microviscosity from ASD simulation results for both CF and CV. 


for both CF and CV conditions shows very little force thickening, even less than the mild 
shear thickening seen in the simulated shear flow. 


8 Summary 

We have used the discrete-particle simulation method Accelerated Stokesian Dynamics (ASD) 
to evaluate the microstructure (pair distribution function, g{Y)) and rheology (viscosity) of 
Brownian colloidal suspensions in active microrheology (MR). In this technique, a probe 
particle is pulled through the suspension bath. The work has considered a wide range of 
flow conditions in terms of the Peclet number Pe, and different volume fractions. This is, to 
our knowledge, the hrst simulation work that considers detailed hydrodynamic interactions 
in active MR. We have considered pulling a probe particle with hxed external force (CF) and 
an approximation to hxed velocity (CV) as analogs to optical tweezer and magnetic bead 
experiments, respectively. Through simulation we demonstrate that the major difference 
between these two conditions at Pe 1 results from the dispersive motion of the probe 
under CF conditions; this dispersion is a result of velocity huctuations associated with the 
heterogeneity of the environment encountered by the probe as well as thermal huctuations, 
with the former dominant at elevated Peclet number for the concentrations studied here. By 
contrast, moving the probe with a hxed velocity, the CV condition, results in a straight-line 
trajectory. As a consequence, the pair distribution function in CF demonstrates one major 
difference compared to CV: when pulling with constant velocity, the length of the zone de¬ 
pleted of bath particles, or wake, formed behind the probe at Pe 2> 1 continuously increases 
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with Peu to the point that it spanned the entire simulation box; see hgure 10. In contrast, 
the ASD simulations for CF conditions show saturation of the length of the wake (and other 
far-held structure) with Pef. We demonstrate that the force-induced dispersion is a direct 
result of many-body hydrodynamic interactions (HI). This effect has not been considered 
in previous theories of active microrheology and simulations since they either ignore HI or 
assume ^ 1 and limit interactions to two-particle level. 

We corroborate the ASD simulation results with a microscopic theory based on the exten¬ 
sion of our theoretical framework for the microstructure and rheology of sheared suspensions 
[Nazockdast & Morris, 20126, 2013]. The theory computes ^'(r) as a solution to the pair 
Smoluchowski equation and is applicable to suspensions with (j) < 0.40 over the entire range 
of Pe. The predicted microstructure g{r) has been used here to compute the average mobility 
of the probe and from this the apparent microrheologically-determined viscosity as a func¬ 
tion of (j) and Pe. A novel aspect of the theory involves modeling the effect of force-induced 
dispersion noted just above to be a factor in CF conditions. For this, we followed the sim¬ 
ple practice adopted in our work on sheared suspensions [Nazockdast &: Morris, 20126] and 
modeled the force-induced relative dispersion flux as a diffusion in the pair Smoluchowski 
equation. 

We were not able to obtain statistically meaningful measure of structural distortion in 
near equilibrium limit, Pe 1 from ASD simulations. Thus we have only used our theory 
to study structure and viscosity in this limit. We hnd that the distorted microstructures 
are quite similar for both hxed force and velocity conditions. This is a consequence of 
the many-body interactions which result in equilibrium structure and a potential of mean 
force, and thus is not in agreement with the predictions of dilute theories, which predict a 
distortion twice as large for constant velocity compared to constant force. We demonstrate 
that including many-body interparticle interactions results in a signihcantly sharper decay of 
structural distortion away from the probe, compared to the predictions of the dilute theory 
at volume fractions as low as (j) = 0.20. As a result of this, the structure and viscosity 
are primarily set by near-contact dynamics of the probe with the bath particles. We show 
that these interactions are very similar in both CF and CV modes of moving the probe, 
resulting in qualitatively similar structure and viscosity. We also hnd that the contribution 
of Brownian forces to the total viscosity in active microrheology is markedly smaller, by 
almost an order of magnitude, than in shear how. 

The theoretical predictions of ^'(r) are in good agreement with ASD results, both near¬ 
contact and at larger distances, for the range of Pe and (j) studied here. Most notably, as 
result of including the ehect of force-induced diffusion, the theory can successfully predict 
the saturation of the wake other structural features away from contact at Pef 2> 1. Our 
simulation results and predictions show that near-contact pair microstructure scales similarly 
with Pe and (j) for CF and CV. We show that the structure at Pe ^ 1 is primarily set by 
the force consumed to pull the particle. This observation is explained by noting that the 
force-induced dihusivity becomes negligible for near contact conhgurations; thus the ratio of 
convection to dihusion huxes scales similarly for both CF and CV pullings as Pe~^ with no 
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strong dependency on cj). As a result the values for viscosity at Pe 3> 1 are quantitatively 
close for CF and CV in predictions and ASD results. In both simulation and theory we 
hnd a mild increase of viscosity with the external force/velocity. The extent of this force¬ 
thickening behavior is, however, smaller than shear thickening in shear flow simulated by 
Stokesian Dynamics. 

We have limited our study here to the case of equal size of probe and bath particles. 
Previous studies have shown a strong influence of probe to bath particle size ratio on mi¬ 
crostructure and microrheology [Almog &: Brenner, 1997]. This theory can be extended to 
consider different size ratios through modihcation of the pair resistance and mobility func¬ 
tions. In this study of active microrheology only hard-sphere suspensions were considered 
whereas the theory can be easily extended to colloidal suspensions with other types of inter¬ 
actions. 
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